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ABSTRACT 

We investigate the location of an ultra-hard X-ray selected sample of AGN from the 
Swift Burst Alert Telescope (BAT) catalog with respect to the main sequence (MS) of 
star-forming galaxies using Herschel-hased measurements of the SFR and M*’s from 
Sloan Digital Sky Survey (SDSS) photometry where the AGN contribution has been 
carefully removed. We construct the MS with galaxies from the Herschel Reference 
Survey and Herschel Stripe 82 Survey using the exact same methods to measure the 
SFR and M* as the iSwi/t/BAT AGN. We hnd a large fraction of the Swift/diPGT 
AGN lie below the MS indicating decreased specihc SFR (sSFR) compared to non- 
AGN galaxies. The Swi/t/BAT AGN are then compared to a high-mass galaxy sample 
(COLD GASS), where we find a similarity between the AGN in COLD GASS and the 
S'wz/l/BAT AGN. Both samples of AGN lie firmly between star-forming galaxies on 
the MS and quiescent galaxies far below the MS. However, we hnd no relationship 
between the X-ray luminosity and distance from the MS. While the morphological 
distribution of the BAT AGN is more similar to star-forming galaxies, the sSFR of 
each morphology is more similar to the COLD GASS AGN. The merger fraction in 
the BAT AGN is much higher than the COLD GASS AGN and star-forming galaxies 
and is related to distance from the MS. These results support a model in which bright 
AGN tend to be in high mass star-forming galaxies in the process of quenching which 
eventually starves the supermassive black hole itself. 
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- infrared: galaxies 


1 INTRODUCTION 

The link between supermassive black holes (SMBH) and 
their host galaxies has been evident for many years through 
the study of correlations between large scale host galaxy 
properties and SMBH mass. Tight correlations were found 
with bulge stellar velocity dispersion (Ferrarese & Merritt 
2000; Gebhardt et al. 2000; Giiltekin et al. 2009), bulge 
luminosity (e.g. Magorrian et al. 1998), and bulge mass 
(Kormendy & Richstone 1995; Marconi & Hunt 2003; 
Haring & Rix 2004) suggesting a coevolution of the SMBH 
with the host galaxy. The question that remains however 
is how two seemingly disjoined objects can influence each 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
t Email: tshimizu@astro.umd.edu 


other over an extremely large range of physical scales. Using 
a simple estimate of the sphere of influence for a SMBH 
(csph = Peebles (1972)) and a typical SMBH 

mass of 10® Mq and stellar velocity dispersion, a — 200 
km s“^, Tsph ~ 10 pc whereas the size of the bulge is 
roughly several kpc (e.g. Simard et al. 2011). Therefore any 
influence from the SMBH must be able to extend over 3 
orders of magnitude in physical scale. 

Active galactic nuclei (AGN), the phase where the 
SMBH is vigorously accreting material and growing, are 
thought to supply the necessary energy to influence the 
galaxy on large scales (e.g. Silk & Rees 1998). This leads 
to a feedback cycle wherein the galaxy supplies cold gas 
that ignites the AGN and fuels star formation, and the 
AGN then returns energy and/or momentum to the galaxy 
that shuts off both accretion and star formation. The ex¬ 
plicit feedback mechanism that runs this cycle is cur¬ 
rently not well understood but is thought to be either 
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large scale outflows (Kaviraj et al. 2011; Clmatti et al. 2013; 
Vellleux et al. 2013; Harrison et al. 2014) or radio jets Best 
(2007); Dubois et al. (2013) (for a review see Fabian 2012). 
Theoretical simulations have also shown that adding AGN 
feedback reproduces well the observed mass and luminosity 
functions while an absence produces too many blue, high 
mass galaxies (e.g. Croton et al. 2006). 

Early evidence for the quenching of star formation 
due to AGN came from studying the colors and stel¬ 
lar masses of large samples of both non-AGN and AGN 
host galaxies. Whereas the non-AGN sample clearly sepa¬ 
rates into two populations on a color-magnitude or color- 
mass diagram, one with red colors (i.e. red sequence 
thought to be quiescent galaxies) and one with blue (i.e. 
blue cloud thought to be strongly star-forming galaxies) 
(Strateva et al. 2001), AGN host galaxies were found to 
be concentrated between them, displaying “green” colors 
(e.g Nandra et al. 2007; Silverman et al. 2008; Hickox et al. 

2009) . Initially, this was interpreted as AGN preferentially 
occurring in galaxies that have had their star formation 
recently quenched (Salim et al. 2007; Martin et al. 2007; 
Schawinski et al. 2009) and are in transition from the blue 
cloud to the red sequence. Optical colors, however, can be 
imprecise tracers of star formation, especially in the presence 
of strong intrinsic dust absorption (e.g. Cardamone et al. 

2010) that obscures recent star formation and causes a red¬ 
dening of the colors that is not due to a reduction of star- 
formation. Green colors can also just be an indication of a 
mixture of old and new stellar populations and not neces¬ 
sarily a “transition” between the two population. Further, 
recent studies that mass-match non-AGN galaxy samples to 
AGN host galaxies reveal that the difference in optical col¬ 
ors virtually disappear (Silverman et al. 2009; Pierce et al. 
2010; Rosario et al. 2013a,b) 

Far-infrared (FIR) emission (A > 40p.m) is essentially 
immune to reddening effects while also being a direct tracer 
of recent star formation. Dust in the galaxy is heated by UV 
photons from recently formed OB stars, that then reemit 
in the mid-far infrared regime (Draine 2003) creating a 
strong correlation between the FIR luminosity and SFR of a 
galaxy (Kennicutt 1998; Kennicutt & Evans 2012). Eurther- 
more, AGN are not thought to strongly affect the FIR (e.g. 
Netzer et al. 2007), and thus the FIR is more robust com¬ 
pared to other SFR indicators such as UV continuum and 
Ha line emission that are mainly used in non-AGN galaxy 
studies. 

The Herschel Space Observatory (Pilbratt et al. 2010) 
opened a window into the FIR universe with the un¬ 
precedented sensitivity of both the Photodetector Array 
Camera and Spectrometer (PACS; Poglitsch et al. (2010)) 
and Spectral and Photometric Imaging Receiver (SPIRE; 
Griffin et al. (2010)) instrument extending the broad band 
spectral energy distributions out to 500 pm and allowing 
an accurate estimate of the FIR luminosity of more objects 
than allowed by previous FIR telescopes (i.e. Infrared Astro¬ 
nomical Satellite (IRAS); Neugebauer et al. (1984)). In this 
study we utilize Herschel to measure the SFRs of a large 
and relatively unbiased sample of AGN and compare their 
location on the SFR-stellar mass (M*) diagram with that of 
the general star-forming galaxy population, which forms a 
“main sequence”. 

The main sequence is the observed tight correlation 


between the stellar mass and SFR of a normal star¬ 
forming galaxy and has been confirmed in depth by many 
studies (e.g. Brinchmann et al. 2004; Noeske et al. 2007; 
Elbaz et al. 2007; Rodighiero et al. 2010; Elbaz et al. 2011; 
Whitaker et al. 2012; Magnelli et al. 2014). This correlation 
seems to exist up to at least z 2 (e.g. Elbaz et al. 2011; 
Whitaker et al. 2012) and possibly all the way to z ~ 4 
(Bouwens et al. 2012; Heinis et al. 2014), with only the nor¬ 
malization changing as a function of redshift, shifting to 
higher SFRs at earlier epochs. This discovery has changed 
theories of galaxy evolution from one that is merger-driven 
to one that is driven more by internal secular processes. 
Elbaz et al. (2011) showed that galaxies that live above the 
main sequence are much more compact with a higher SFR 
surface density, indicative of a major merger. Main sequence 
galaxies, though, have disk-like morphologies inconsistent 
with a recent merger (Wuyts et al. 2011) that suggests star 
formation is triggered by internal processes such as disk in¬ 
stabilities. 

In the past, the question of where AGN fit into the pic¬ 
ture was unclear because most main sequence studies pur¬ 
posely excluded AGN due to its messy contribution to SFR 
indicators. However, with Hersehel, more accurate estimates 
of the SFRs for AGN host galaxies can be calculated as well 
as better detection rates. For example, both Mullaney et al. 
(2012b) and Rosario et al. (2013b), using deep observations 
of large fields examined the star-forming properties of X-ray 
selected AGN. Both came to the conclusion that AGN pri¬ 
marily reside in main sequence galaxies, calling into question 
the long-held idea that AGN host galaxies are in the pro¬ 
cess of quenching star formation. If most AGN are in main 
sequence galaxies, this could indicate that moderate lumi¬ 
nosity AGN are simply coincidental with a large cold gas 
reservoir that also fuels star formation (Vito et al. 2014). 

In this Paper, we seek to fully investigate the location 
of AGN on the SFR-M* plane using Herschel observations 
of Swift/BAT selected AGN. These are ultra-hard X-ray 
confirmed AGN at low redshift using a selection method 
that is unbiased with respect to both obscuration and host 
galaxy contamination. We calculate SFRs using Herschel 
photometry and a simple model to disentangle star forma¬ 
tion and AGN contributions and combine them with esti¬ 
mates of the stellar mass using AGN subtracted SDSS pho¬ 
tometry from Koss et al. (2011). The SFRs and M*’s are 
then compared to a local normal star-forming galaxy sample 
to define the main sequence as well as a sample of galaxies 
purely selected on stellar mass. Finally we discuss the impli¬ 
cations of our results and compare them to previous stud¬ 
ies. Throughout this paper we use a ACDM cosmology with 
Ho = 70kms“^ Mpc^, = 0.3, and Ha = 0.7. Luminosity 
distances were calculated using this cosmology along with 
redshifts taken mainly from the NASA/IPAC Extragalactic 
Database (NED)^, except for those objects with z < 0.01 
where we used measured distances from the Extragalactic 
Distance Database^. We either use or correct for a Chabrier 
or Kroupa initial mass function (IMF) for all star formation 
rate calculations. 


^ http://ned.ipac.caltech.6du/ 
^ http://6dd.ifa.hawaii.edu/ 
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2 SAMPLES AND OBSERVATIONS 


2.1 Swift/BAT AGN 

Our parent sample of AGN was drawn from the 58 month 
Swift Burst Alert Telescope (BAT) (Gehrels et al. 2004; 
Barthelmy et al. 2005) catalog (Baumgartner et al. 2013) 
with a redshift cutoff of z < 0.05, totaling 313 AGN (149 
Seyfert 1-1.5s, 157 Seyfert 1.8-2s, 6 LINERs, and 1 uniden¬ 
tified AGN). The catalog is the result of continuous mon¬ 
itoring by Swift/BAT of the entire sky in the 14“195 keV 
energy range. These high energies allow for an unambiguous 
detection of AGN with little to no contamination from the 
host galaxy and significantly reduced selection effects due to 
obscuration. 

All 313 AGN were observed by the Herschel Space Ob¬ 
servatory with 291 part of our program (PI: R. Mushotzky, 
PID: OTl_rmushotz_l) and the remaining 22 obtained 
from other programs publicly available on the Herschel Sci¬ 
ence Archive. The sample was imaged by both the PACS 
and SPIRE instruments providing, for the first time, sensi¬ 
tive FIR photometry from 70-500 pm for a large, ultra-hard 
X-ray selected sample of AGN. Detailed descriptions of the 
reduction and analysis of the PAGS and SPIRE images are 
given in Melendez et al. (2014) and Shimizu et al (2015, in 
preparation), but we provide a brief summary here. 

PACS and SPIRE together imaged the sample in 5 
broad bands: 70 and 160 pm (PACS) and 250, 350, and 500 
pm (SPIRE). Level 0 (raw) data were reduced to Level 1 
using the standard pipeline provided by the Herschel In¬ 
teractive Processing Environment (Ott 2010) v8.0. Maps 
were produced from the Level 1 data using Scanamorphos 
(Roussel 2013) vl9.0, a software package that takes advan¬ 
tage of the redundancy inherent in the scanning procedure 
of Herschel to remove both thermal and non-thermal low- 
frequency noise. Circular and elliptical apertures with radii 
chosen visually to encompass the entirety of the FIR emis¬ 
sion were then used to extract the photometry for each wave¬ 
band. la errors for the photometry were determined using 
a combination of the pixel-by-pixel errors in the aperture, 
an estimate of the root-mean-square of the background, and 
calibration uncertainty. 

In addition to Herschel observations, Koss et al. (2011) 
analyzed optical images of 185 BAT AGN from the Sloan 
Digital Sky Survey (SDSS) and Kitt Peak National Obser¬ 
vatory. Using GALFIT (Peng et al. 2002), they were able 
to accurately measure the host galaxy light by subtract¬ 
ing out the central point source due to the AGN. Reliable 
stellar masses for the BAT AGN host galaxies were then es¬ 
timated using standard stellar population models. Because 
the Koss et al. (2011) BAT sample was chosen from the 22 
month catalog (Tueller et al. 2010), 45/185 were not ob¬ 
served with Herschel, reducing the sample to 140 AGN. 
Furthermore, Koss et al. (2011) flagged 18/140 objects for 
incomplete PSF-subtraction from the griz images, so we 
choose not to include these sources resulting in a final sample 
of 122 AGN, including 46 Sy Is, 72 Sy 2s, and 4 LINERs, 
where we define a Sy 1 as Sy 1-1.5 and Sy 2 as Sy 1.8-2. 
The reason for the discrepancy between the number of Sy 
Is and Sy 2s is that all 18 of the objects that were flagged 
for incomplete PSF-subtraction are Sy Is. 


2.2 Herschel Reference Survey 

To form the main sequence, we need a large and complete 
sample of star-forming galaxies that do not host an AGN 
but have been observed at the same wavelengths allowing 
for a consistent determination of both the SFR and stel¬ 
lar mass. For these reasons we chose the Herschel Reference 
Survey (HRS; Boselli et al. (2010)), a guaranteed time Her¬ 
schel key project that imaged 323 K-band selected galaxies 
from 100-500 pm. The HRS spans all morphological types 
and was volume limited to contain galaxies between 15 and 
25 Mpc away. Even though our sample stretches out to 
z = 0.05 (~200 Mpc), the HRS represents the best sample 
of local star-forming galaxies to compare with given that 
both have been observed by Herschel as well as other tele¬ 
scopes including the Wide-Field Infrared Survey Explorer 
(Wright et al. 2010, WISE), the Galaxy Evolution Explorer 
(GALEX Martin et al. 2005), and SDSS. 

The HRS PACS and SPIRE images were analyzed in 
Cortese et al. (2014) and Ciesla et al. (2012) producing pho¬ 
tometry at 100, 160, 250, 350, and 500 pm. We applied the 
same corrections to the SPIRE photometry as described in 
Cortese et al. (2014) due to changes in the SPIRE calibra¬ 
tion and beam size. The corrections reduce the SPIRE flux 
densities by 7, 6, and 9 per cent at 250, 350, and 500 pm 
respectively. 

A subset of the HRS are galaxies within the Virgo Clus¬ 
ter and have been affected by the dense environment through 
the stripping of their gas (Boselli et al. 2006). The dust con¬ 
tent of these galaxies has also been shown to be affected by 
the environment (Cortese et al. 2010, 2012a). Therefore, fol¬ 
lowing Ciesla et al. (2014), we restrict the HRS sample to 
only those galaxies that aren’t “Hl-deficient” as defined in 
Boselli et al. (2012) which reduces the sample to 146 galax¬ 
ies. 


2.3 COLD GASS 

While the HRS represents a sample that was observed using 
the same telescopes, it is limited in its range of stellar mass 
especially above 10^° Mq. For a complete comparison to 
the BAT AGN, we supplemented HRS with the CO Legacy 
Database for GASS (COLD GASS; Saintonge et al. 2011a), 
a 366 galaxy subsample of the GALEX Arecibo SDSS Survey 
(GASS; Catinella et al. 2010). The GASS sample consists of 
~1000 galaxies randomly selected such that every galaxy lies 
within the footprint of the SDSS spectroscopic survey, the 
ALFALFA survey, and the GALEX Medium Imaging Sur¬ 
vey. The galaxies also were selected to have a redshift range 
0.025 < z < 0.05 and a stellar mass between 10^° and 10^^'®, 
both of which match very well to the BAT AGN. 366 galax¬ 
ies were randomly selected from GASS to form the COLD 
GASS sample and have deep CO(l-O) imaging with IRAM. 
COLD GASS then represents a completely unbiased sample 
of galaxies above 10^° Mq. 

As described in Catinella et al. (2010) and 
Saintonge et al. (2011a), both GASS and COLD GASS 
were selected to have a uniform logM* distribution. The 
stellar mass distribution that is observed, however, is more 
heavily weighted towards lower mass galaxies (see Figure 
1 of Saintonge et al. (2011a)). Therefore, Saintonge et al. 
(2011a) constructed 50 representative subsamples of COLD 
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GASS that matches the observed M* distribution. Each 
subsample contains between 200-260 galaxies and can be 
used to test the robustness of any relation that might 
be observed. Throughout this Paper, we explicitly note 
whether the full COLD GASS sample, an average of 
the representative subsamples, or a single representative 
subsample is being used. 

Because COLD GASS was selected in an unbiased way, 
the sample contains a mixture of galaxy types from star¬ 
forming to quiescent to AGN. To determine the type of 
each galaxy in COLD GASS, we cross-matched the sam¬ 
ple with the MPA-JHU SDSS DR7 spectroscopic catalog.® 
In this catalog, galaxies that have SDSS optical spectra were 
classified according to their location on the standard BPT 
diagram (Baldwin et al. 1981). Brinchmann et al. (2004) de¬ 
scribes the details of the classification^. Galaxies were sep¬ 
arated into 6 groups: Star-forming, Low S/N Star-forming, 
AGN, Composite, LINER, and Quiescent. For the purposes 
of comparison with the BAT AGN, we combined the Star¬ 
forming and Low S/N Star-forming groups into a single 
Star-forming group and the AGN and Composite groups 
into a single AGN+Composite group. We combine the AGN 
and Composite group because after cross-matching the BAT 
AGN with the same SDSS DR7 sample, we find that al¬ 
most all are classified as an AGN or Composite galaxy. The 
COLD GASS AGN+Composite subsample then represents 
an ideal optically selected sample of AGN to compare with 
the BAT AGN, while the Star-forming and Quiescent sam¬ 
ple represent ideal non-AGN samples. We denote the differ¬ 
ent subsamples In the following way: Star-forming as CCS, 
AGN+Composite as CCA, LINERs as CGL, and Quiescent 
as CGQ. 


2.4 Herschel Stripe 82 Survey 

With the HRS and COLD GASS samples we have one that 
matches the photometry available to the BAT AGN (HRS) 
and one that more closely matches the physical properties 
of the BAT AGN host galaxies (COLD GASS). HRS lacks 
galaxies at the high stellar mass of the BAT AGN whereas 
COLD GASS lacks Herschel and WISE photometry to allow 
for a consistent comparison between SFRs. Therefore, we use 
a third comparison sample to the BAT AGN, the Herschel 
Stripe 82 survey (HerS; Viero et al. 2014)that bridges the 
gap between HRS and COLD GASS. 

HerS is a 79 deg^ survey of SDSS Stripe 82 with Her¬ 
schel/SPIRE at 250, 350, and 500 qm. Viero et al. (2014) 
produced the catalog of HerS sources while Rosario et al 
(2015, in preparation) calculated SFRs based on FIR lumi¬ 
nosities. Given its small area, the volume covered at low 
redshift is smaller than that of both the COLDGASS and 
BAT samples which means HerS is rather incomplete at high 
stellar masses (> 10.5Mq). However it is still better than 


® http://www.mpa-garching.mpg.de/SDSS/DR7/ and http:// 
home. strw. leideniiniv.nl/~jarle/SDSS/ 

We note that Saintonge et al. (2012) used a more conservative 
method to classify AGN in the COLD GASS. A galaxy was con¬ 
sidered an AGN if log[NII/HQ]> —0.22 and log[OIII/H/3]> 0.48. 
This led to a much lower fraction of AGN in their analysis (6 per 
cent) compared to this work (30 per cent). 


the HRS and has Herschel photometry available, albeit only 
from the SPIRE instrument. 

We selected our HerS sample from Rosario et al (2015, 
in preparation), who matched HerS sources from Viero et al. 
(2014) with the SDSS MPA-JHU DR7 catalog. Because the 
HerS catalog assumes all galaxies are point sources, Rosario 
et al (2015, in preparation) limited their sample to 2 > 0.02. 
We further limit the sample to 2 < 0.08 to match the BAT 
AGN while also pushing out to a slightly larger volume to 
populate the high mass end better without a significantly 
affecting our results. After cutting sources which do not have 
a measured stellar mass as well as sources with low S/N 
(< 3) emission lines (only 5 objects) this results in a final 
HerS sample of 517 objects. Due to the combined magnitude 
cut from the SDSS and the limited sensitivity of SPIRE from 
which the HerS catalog was built, Malmquist bias can be 
a problem. However, because of the relatively low redshift 
nature of our sample we do not expect it to largely bias our 
results. 

The HerS sample was also split into some of the same 
classifications as COLDGASS using the BPT diagram. HerS 
contains both a star-forming and AGN population and will 
be designated as HerS SF and HerS AGN. The HerS AGN, 
just as CGA, is a combination of the AGN and Compos¬ 
ite classifications. Within HerS, there is also an Uncertain 
classification which indicates a galaxy that is missing or has 
upper limits for at least one of the 4 key lines needed to 
classify it using the BPT diagram. 


3 STAR FORMATION RATE ESTIMATION 

The unique Herschel data provides a means for determining 
accurate star formation rates (SFR). IR emission has long 
been used as a calibrator for star formation (see Kennicutt 
1998, for a review), because it probes the dust population 
that reprocesses the UV emission from young stars (e.g. 
Draine 2003). The specific wavelength range of Herschel cov¬ 
ers the bulk of the IR emission from dust including the char¬ 
acteristic FIR bump typically seen in star-forming galaxies 
(e.g. Dale et al. 2007) allowing for precise measures of the 
total IR luminosity and thus the star formation rate, espe¬ 
cially for AGN host galaxies where many of the often used 
SFR indicators (e.g. Ho, UV continuum) can be substan¬ 
tially contaminated by AGN-related emission. 

Many SED-fitting packages exist in the literature rang¬ 
ing from template based models to full dust radiative trans¬ 
fer models. However, given the low number of data points for 
our SEDs (at most seven), we chose to fit our SEDs with the 
model described in Casey (2012), which is a combination of 
an exponentially cutoff mid-infrared (MIR) power law and 
a single temperature greybody. Details and results of the 
SED fitting for the BAT AGN will be given in a forthcom¬ 
ing paper (Shimizu et al 2015, in preparation), however we 
provide example fits to both the BAT AGN and HRS in Ap¬ 
pendix A as well as a brief overview of the fitting procedure. 
In Section 6.1 we discuss extensively and test whether this 
model introduces systematic biases especially related to the 
decomposition of the SED. 

The Herschel photometry constrains the greybody com¬ 
ponent, but we need additional shorter wavelength data to 
constrain the power law component. Therefore, we cross- 
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correlated our sample with the Wide-field Infrared Survey 
Explorer (Wright et al. 2010, WISE) All-Sky Release Cata¬ 
log on the IRSA website.® We only use the 12 (W3) and 22 
(W4) pm photometry to avoid contamination by the stel¬ 
lar population in the host gala:xy. WISE photometry for the 
HRS were taken from Ciesla et al. (2014). 

For consistency, we fit both the BAT AGN and HRS 
gala:xies using the Casey (2012) model even though the HRS 
gala:xies do not host an AGN or are classified as ULIRGs. In 
this way, we can account for a portion of MIR emission that 
is due to normal star formation rather than AGN heated 
dust. We fit all BAT AGN and HRS galaxies with at least 
four detected photometric points. This restriction ensures 
quality photometry for each galaxy and removes only a fur¬ 
ther 9 and 11 galaxies from the BAT AGN and HRS sam¬ 
ple respectively for final sample sizes of 113 and 135. We 
determine three luminosities for each galaxy: a total IR lu¬ 
minosity (I/tir), a MIR power law luminosity (Lmir) and 
a greybody luminosity (Lcrey). Each one was calculated by 
integrating the best-fitting model from 8-1000 pm. Ltir is 
the luminosity from integrating over the total model while 
Tmir is from only integrating the MIR power law compo¬ 
nent, and Lorey is from the greybody component. 

Star formation rates are then calculated using one of 
these IR luminosities and the equation from Murphy et al. 
( 2011 ). 


SFRir = 


Lir [erg s 
2.57 X 10« 


( 1 ) 


For the HRS galaxies Lir = Ltir since there is no AGN to 
contribute to the IR emission. For the BAT AGN, however, 
we use Lir = 4/3LGrey. The 4/3 is a correction factor to 
account for MIR emission from star formation. To determine 
it, we calculated the average ratio of LMm/Lcrey for the 
HRS sample, which contain no AGN, and found it to be 
narrowly distributed around 1/3. This means that only using 
Lcrey to determine the SFR underestimates it by 1 /3 so we 
need to multiply Lcrey by 4/3 as a correction. 

While FIR emission probes dust obscured star forma¬ 
tion, the UV continuum is a measure of unobscured star for¬ 
mation by tracing the direct light from young massive stars. 
Hence, a complete census of star formation can be found by 
combining measurements from both wavebands. AGN are 
strong emitters in the UV though, so only Seyfert 2 galax¬ 
ies will have reliable UV photometry to combine with the 
FIR for a SFR. Using GALEX far-UV data from the GCAT 
catalog (Bianchi et al. 2014) for the BAT AGN Seyfert 2’s 
and the same from Boselli et al. (2013) for the HRS, we cal¬ 
culated dust-corrected UV SFRs. We found that using UV 
SFRs for both the HRS and BAT AGN Seyfert 2s had no ef¬ 
fect on the results of this Paper. Thus, we choose to use the 
FIR only SFRs to allow a larger BAT AGN sample (Seyfert 
I’s and 2’s). 

The COLD GASS sample unfortunately was not ob¬ 
served with Herschel and does not allow for the same 
calculation of the SFR. We use the SFRs provided in 
Saintonge et al. (2011b) which were calculated by fit¬ 
ting SDSS and GALEX photometry to Bruzual & Chariot 
(2003) models. Saintonge et al. (2011b) compared these 


® http://irsa.ipac.caltech.edu/Missions/wise.html 


SFRs to those inferred from combined GALEX FUV and 
Spitzer 70 pm data finding a strong correlation with only 
a scatter of 0.22 dex and indicating FIR-based SFRs are 
consistent with optical-UV ones. We recognize that the 
Saintonge et al. (2011b) comparison however did not use any 
Hersehel photometry, but three calibration analyses®^® show 
that Spitzer and Hersehel produce consistent fluxes. Further 
Dominguez Sanchez et al. (2014) recently showed that SFRs 
calculated from FIR SED fitting are consistent (with large 
scatter especially for pure AGN and quiescent galaxies) with 
the SFRs inferred from the MPA-JHU SDSS spectral anal¬ 
ysis. Based on all these indirect tests, we are confident that 
the COLD GASS SFRs are consistent with the Herschel 
based ones we use for the HRS and BAT AGN in a way that 
does not effect the results of this Paper. 

For the HerS sample, the SFRs were estimated by fit¬ 
ting the SPIRE 250 pm and WISE 22 pm fluxes to the 
Dale & Helou (2002) (DH02) templates. Each of the 64 tem¬ 
plates represents a different value of a where a is the power 
law index for the distribution of dust mass over heating in¬ 
tensity in a galaxy, minimization was used to scale each of 
the 64 templates to the observed fluxes in the HerS sample, 
then the template with the lowest was chosen as the best 
fit. The best-fit template was integrated between 8-1000 pm 
to calculate Lir and converted to a SER using the SER-L/r 
relation from Kennicutt (1998) adjusted for a Chabrier IMF 
(lowered by a factor of 1.7). In the absence of WISE 22 pm 
photometry, a was fixed at -2.0. 

This is a slightly different method than the one used 
for the HRS and BAT AGN. To test for systematics we ht 
the HRS galaxies with 250 pm and 22 pm detections using 
the DH02 method and compared the SFRs. A linear fit to 
the two SFRs reveals a slope of 1.0 and an offset of 0.11 
dex with the SFRs determined from the Casey (2012) and 
Equation 1 higher than the ones from DH02 and the relation 
from Kennicutt (1998). Therefore we adjust all of the HerS 
SFRs by adding a constant of 0.11 dex. 


4 STELLAR MASS ESTIMATES 

SFRs are only one-half of the main sequence; stellar masses 
are also needed for the galaxies. Cortese et al. (2012b) cal¬ 
culated the stellar masses of the HRS using the relation from 
Zibetti et al. (2009): 

log(M*/Li) =-0.963 -bl.032(g-f) (2) 

where M* is the stellar mass and Li is the i-band lumi¬ 
nosity, both in solar units. To be consistent we also used 
this equation for the BAT AGN. The g — i color was calcu¬ 
lated using the PSF-subtracted photometry from Koss et al. 
(2011). The stellar masses determined here correlate very 
well with the stellar masses from Koss et al. (2011) with 
a Pearson correlation coefficient, rp = 0.85 indicating a 


® http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/ExtSrcPhotom.pdf 
^ http://herschel.esac.esa.int/twiki/bin/ 
viewfile/Public/PacsCalibrationWeb?rev=l; 
filenaiiie=PlCC-NHSC-TN-029. pdf 

® https://nhscsci.ipac.caltech.edu/pacs/docs/Photometer/ 
PlCC- NHSC- TR- 034. pdf 
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highly linear relationship. However the Koss et al. (2011) 
values are systematically larger by a factor of ~ 2® Since 
the goal of this Paper is to compare AGN with normal 
star-forming galaxies and not absolute measures of stel¬ 
lar mass and SFR, we choose to use Equation 2 for the 
BAT AGN stellar masses. We performed the same calcula¬ 
tion for the GOLD GASS galaxies and compared these M* 
with those provided in Saintonge et al. (2011b). We find the 
same strong correlation and the same systematic offset as the 
BAT AGN stellar masses so we choose to also use the stellar 
masses calculated in this Paper for GOLD GASS as well. The 
HerS stellar masses are taken from the MPA-JHU database 
which used the same method as both Koss et al. (2011) and 
Saintonge et al. (2011b). However, we do not have g ox i 
photometry for HerS so we apply a factor of 2 correction to 
them to match the stellar mass scale from Equation 2. 

5 RESULTS 

5.1 Location of AGN in SFR-M* Plane 

One key issue in this analysis is which main sequence to use. 
Many authors have published main sequence relations (for a 
nice compilation see: Speagle et al. 2014), however each re¬ 
lation was determined differently using different stellar mass 
estimates, SFR indicators, and redshift ranges. This has re¬ 
sulted in a large spread of values for both the slope and 
normalization of the main sequence, especially in the local 
universe. Therefore, we choose to calculate our own main 
sequence relation using only the HRS galaxies and the HerS 
star-forming galaxies since both the stellar masses and SFRs 
were calculated with comparable methods as for the BAT 
AGN. We use a linear bisector (Isobe et al. 1990) to fit the 
HRS+HerS data, resulting in the following equation: 

log(SFRiR) = 1.01 log(M.) - 9.87 (3) 

The scatter of the relation is 0.36 dex, similar to the scat¬ 
ter seen in other studies (e.g Noeske et al. 2007; Peng et al. 
2010). The slope and the normalization are slightly steeper 
and smaller respectively than that found in Peng et al. 
(2010), who analyzed the SFRs and stellar masses of the 
entire local (0.02 < 2 < 0.085) SDSS population. The slope 
is much steeper than the 2 = 0 slope predicted using the 
Speagle et al. (2014) relation (0.5). This is possibly due to 
the addition of lower mass objects from HRS as well differ¬ 
ences in the measurement of the SFR and M*. Further, due 
to the large scatter in MS relations measured at low red- 
shift, Speagle et al. (2014) specifically did not include low 
redshift studies in formulating their redshft-dependent MS 
relation. This is the reason we set out to formulate our own 
main sequence relation that uses a well-defined star-forming 
galaxy sample and measures of the SFR and stellar mass 
that are consistent between the non-AGN and AGN host 
galaxy samples. 


® This is most likely due to the different stellar population models 
used in Zibetti et al. (2009). They used the 2007 version of the 
Bruzual & Chariot (2003) models that included a new presciption 
for thermally pulsing AGB stars. This decreases the stellar mass 
by a factor of 2 especially for star-forming galaxies. This has no 



A log SFR 

Figure 2. Kernel density estimate (KDE) of the logarithmic dis¬ 
tance (A log SFR) from the main sequence for the BAT AGN, a 
representative subsample of COLD GASS, and the HerS sample. 


Figure 1 plots the HRS, BAT AGN, a randomly cho¬ 
sen representative subsample of COLD GASS, and the HerS 
sample on the SFR-M* plane along with Equation 3 and 
its scatter. Visually it is clear that a large fraction of the 
BAT AGN and COLD GASS lie either inside or below the 
MS. Table 1 quantifies the exact fraction of galaxies in five 
different regions. The different regions are divided accord¬ 
ing to A log SFR = log SFRobs — log SFRms where SFRobs 
is SFRir and SFRms is the SFR expected given the M* 
of the galaxy using Equation 3. A log SFR represents the 
distance a source is from the main sequence and given the 
nearly linear MS relation is proportional to specific SFR 
(sSFR = SFR/M*). The five regions are defined as: above 
the main sequence (A log SFR > Icr), inside the main se¬ 
quence (Icr > A log SFR > —Icr), 1-2(7 below the main se¬ 
quence (—1(7 > A log SFR > —2(7), 2-3(7 below the main se¬ 
quence (—2(7 > A log SFR > —3(7), and 3a below the main 
sequence (A log SFR < —3a) with a equal to the observed 
scatter in the main sequence relationships (i.e. 0.36 dex). We 
break the “below” region in three separate regions to judge 
how the sample is clustered. If all of the sources below the 
main sequence are in the Icr region, then it could be argued 
that most of the AGN are main sequence galaxies and sim¬ 
ply display a larger scatter. The fractions for the different 
GOLD GASS subsamples are the average fraction over all 
50 representative subsamples. 

The numbers confirm the visual impression seen in Fig¬ 
ure 1 that the BAT AGN mainly live inside the main se¬ 
quence or below it. 28 per cent lie inside the main sequence 
and 66 per cent lie below it from adding together the three 
“below” regions. Only 5 per cent of the sample is above the 
main sequence. The sources below the main sequence are 
also well spread out between the Icr, 2a, and 3a regions, 
showing the BAT AGN do not cluster near the edge of the 


effect on our results as long as all samples are on the same stellar 
mass scale. 
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M* [Mq] 


Figure 1. The relationship between SFR and M* for the HRS (black dots), BAT AGN (top: colored stars), a representative subsample 
of COLD GASS (middle: colored symbols), and the HerS sample (bottom: colored symbols). The solid line represents the main sequence 
relationship calculated using the IR-based HRS+HerS SFRs (Equation 3). The dashed lines are the measured Icr (0.32 dex) scatter in 
the relationship. 
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Table 1. Main Sequence Locations for BAT AGN and COLD GASS 


Sample 

(1) 

Total 

(2) 

Above MS 
(3) 

Inside MS 
(4) 

1 — 2(j Below 
(5) 

2 — 3(T Below 
(6) 

> 3cr Below 

(7) 

BAT AGN 

All 


0.05 

0.28 

0.28 

0.18 

0.20 

Sy 1 

0.35 

0.05 

0.20 

0.35 

0.20 

0.20 

Sy 2 

0.61 

0.06 

0.30 

0.26 

0.17 

0.20 

LINER 

0.04 

0.00 

0.75 

0.00 

0.00 

0.25 

COLD GASS (CG) 

All 


0.02 

0.33 

0.15 

0.12 

0.39 

Star-Forming (CGS) 

0.37 

0.04 

0.62 

0.12 

0.05 

0.16 

AGN+Comps (CGA) 

0.30 

0.02 

0.28 

0.23 

0.15 

0.34 

LINER (CGL) 

0.17 

0.00 

0.07 

0.18 

0.25 

0.50 

Quiescent 

0.15 

0.00 

0.00 

0.00 

0.03 

0.97 

Herschel Stripe 82 (HerS) 

All 


0.12 

0.67 

0.14 

0.05 

0.02 

Star-Forming 

0.72 

0.14 

0.78 

0.07 

0.01 

0.00 

AGN+Comps 

0.19 

0.08 

0.48 

0.25 

0.10 

0.09 

Uncertain 

0.08 

0.02 

0.21 

0.42 

0.28 

0.07 


Notes (1) Sample and subsample names for the BAT AGN and COLD GASS. (2) Fraction of the total sample (i.e. 
BAT AGN and COLD GASS) that each subsample occupies. (3) Fraction of sample that is above the main sequence 
(MS) (AlogSFR > Icr; O' = 0.32 dex). (4) Fraction of sample that is inside the MS (Icr > AlogSFR > —Icr). (5) 
Fraction of sample that is between Icr and 2cr below the MS (—Icr > AlogSFR > —2cr). (6) Fraction of sample that 
is between 2cr and 3cr below the MS (—2cr > AlogSFR > —3cr). (7) Fraction of sample that is greater than 3cr below 
the MS (AlogSFR < —3cr). The COLD GASS fractions represent the average fractions over all 50 representative 
subsamples. Due to round-off errors, the sums across each row are not exactly equal to 1. 


main sequence and in fact a significant percentage (20 per 
cent) display SFRs more than 3cr below what is expected 
given their stellar masses. 

For the COLD GASS sample as a whole, there seems 
to be a bimodal distribution with 33 per cent of galaxies in¬ 
side the main sequence and 39 per cent > 3a below it. This 
mirrors what has been extensively seen using optical colors 
and is another representation of the split into the blue cloud 
and red sequence. The BAT AGN do not display the same 
bimodality and overall have a very different distribution of 
A log SFR as shown in Figure We ran K-S tests to com¬ 
pare the BAT AGN with each COLD GASS representative 
subsample (see Section 2.3 for a description of the subsam¬ 
ples) and found that 0/50 tests returned a probability > 5 
per cent verifying that the two samples are not drawn from 
the same parent distribution. 

The HerS sample differs from the COLD GASS 
and BAT AGN samples and is sharply peaked around 
A log SFR = 0, signifying most of the galaxies are on the 
MS. Indeed, 67% of the total HerS sample lie inside the main 
sequence, while only 12 per cent are above and 21 per cent 

We take this opportunity to explain that we choose to rep¬ 
resent distributions of values using a Kernel Density Estimate 
(KDE) rather than a histogram. This is due to the visualization 
of a histogram being highly dependent on the bin size, number 
of bins chosen, and the edges of the bins. A KDE represents each 
point in a data set with a specific kernel and sums all of them to¬ 
gether. In this Paper we use a Gaussian kernel. The only tunable 
parameter is the kernel width for which we use “Scott’s Rule” 
(Scott 1992), width = where N is the number of data 

points. 


below. HerS likely does not reach as far below the main se¬ 
quence as GOLD GASS given the parent sample was selected 
based on a SPIRE 250 qm detection. Using the Viero et al. 
(2014) 3a depth for SPIRE 250 qm of 65 mjy and assuming 
a greybody with Tdust = 25 K and p = 2.0, the minimum 
SFR detectable is 0.5 Mq yr“^, much higher than the level 
reached by COLD GASS and the BAT AGN. Even so, we 
again ran a K-S test between the BAT AGN and HerS and 
found a Pk-s « 0.01 indicating they are drawn from dif¬ 
ferent populations. 

It is only when we separate the GOLD GASS and HerS 
sample into their different classifications do we find a simi¬ 
larity. Figure 3 shows the KDEs of A log SFR for the BAT 
AGN and each COLD GASS and HerS classification. The 
COLD GASS KDEs were calculated from a randomly chosen 
representative subsample. The BAT AGN are most similar 
to the CGA (e.g. the AGN population of GOLD GASS), dis¬ 
playing lower values of AlogSFR than CGS/HerS SF and 
higher values than GGL and especially CGQ. In fact the 
percentages in each region for the BAT AGN and CGA are 
nearly identical except in the > 3 ct region where there is 
a larger fraction of GGA. We again ran K-S tests for each 
of the 50 representative subsamples for GOLD GASS and 
found 48/50 CGA subsamples returned a probability > 5 
per cent indicating that the BAT AGN and GGA are con¬ 
sistent with the same parent population. 0/50 of the CGS, 
CGL, and CGQ subsamples returned a probability > 5 
per cent of being consistent with the BAT AGN. The K- 
S test between the BAT AGN and HerS AGN returned a 
Pk-s = 0.002 while the ones with HerS SF and HerS Un¬ 
certain returned Pk-s « 0.001 and Pk-s = 0.07. Using 
the standard 5 per cent cutoff to determine if the samples 
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Figure 3. Same as Figure 2 but with the COLD GASS subsample split into 4 classifications based on their location in the BPT diagram. 
CGS; Star-forming galaxies. CGA: AGN and composite galaxies. GGL: LINERs. CGQ: Quiescent galaxies. The HerS sample was split 
into 3 classifications (Star-forming, AGN and composite, and Uncertain). The BAT AGN galaxies show a similar distribution in A log 
SFR as the CGA and HerS AGN sample (upper right). 


originate from the same population would indicate that the 
BAT AGN are most similar to the HerS Uncertain galaxy 
population. However as we discuss above, this is most likely 
due to the lower SFR depth reached in HerS. Also the HerS 
Uncertain classification is a more conservative classification 
and includes many galaxies that would have been classified 
using the Brinchmann et al. (2004) system including AGN 
and LINERs. The HerS AGN still show lower SFRs than 
the HerS SF sample with 48 per cent inside the MS com¬ 
pared to 78 per cent and 44 percent below it compared to 
8 per cent. Therefore, both an optically selected sample of 
AGN and an ultra-hard X-ray selected sample of AGN dis¬ 
play the same property: they lie in between a strongly 
star-forming and quiescent group and suggests that 
AGN host galaxies are in transition between the two 
populations. 

We note that while the K-S tests confirm that the BAT 
AGN and GGA are from the same population, both Ta¬ 
ble 1 and Figure 3 show that the BAT AGN contain slightly 


more galaxies with higher SFRs while GGA contains slightly 
more quiescent galaxies. We hypothesize this is due to the 
selection method for the two groups. The BAT AGN are 
X-ray selected, a method that is completely independent 
from the star-forming properties of the host galaxy while 
the CGA are selected by optical emission line ratios whose 
origin can be a mixture of AGN and star-formation. If a 
galaxy is highly star-forming, optical emission line ratios 
are more likely to classify it as a star-forming galaxy rather 
than AGN or even composite (see Trump et al. (2015) for 
biases associated with line ratio selection of AGN). Indeed, 
many of the BAT AGN above the main sequence are either 
involved in a merger or are known starburst galaxies (e.g. 
Mrk 18, NGC 3079, NGC 7679). 

Because the BAT survey is flux-limited, it is biased 
against weak AGN, while COLD GASS, due to its selec¬ 
tion from SDSS, can reach to lower AGN luminosities. If 
low-luminosity AGN are more associated with quiescent, 
early-type galaxies (e.g. Kauffmann et al. 2003), this would 
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explain the larger fraction of the CGA group in the > Scr 
region compared to the BAT AGN. In Section 5.4 we show 
that the GGA contain a larger fraction of elliptical galaxies 
than the BAT sample which have lower values of sSFR. 

We also note that the BAT AGN and CGS have essen¬ 
tially the same fraction of galaxies in the > 3 (t below region 
(16 per cent vs. 14 per cent). This is seen in Figure 3 as the 
long tail towards low A log SFR for GGS. However, this does 
not change the general result that the BAT AGN in general 
show lower levels of sSFR than the star-forming sample due 
to the much higher percentages in the 1 — 2g and 2 — 3a 
below regions. BAT AGN occur in the 1 — 2a region at a 
> 3 times higher rate than the CGS and more than 2 times 
in the 2 — 3cr below region. Further 64% of the CGS occur 
inside the main sequence compared to only 40% of the BAT 
AGN. So over the entire population, AGN are more likely 
to be found in host galaxies that have lower SFRs than the 
main sequence. 

5.2 Differences Between Seyfert Is and 2s 

According to the unified model (Antonucci 1993; 
Urry & Padovani 1995), orientation distorts our view 
of AGN and causes the differences seen between Sy Is and 
2s. Current models invoke an anisotropic dusty and possibly 
clumpy torus (e.g. Nenkova et al. 2008; Honig & Kishimoto 
2010) that absorbs and scatters the nuclear optical/UV/X- 
ray emission. Two regions of ionized gas produce optical 
emission lines, one at relatively close distances to the central 
AGN that produces broad emission lines (i.e. broad line 
region, BLR) and one at further distances that produces 
narrow emission lines (i.e. narrow line region, NLR). Under 
the unified model, Sy 1 galaxies, which display very bright 
nuclear point sources as wells as broad optical emission 
lines are viewed along lines of sight through the opening 
angle of the torus allowing access to the BLR and accretion 
disk. Sy 2 galaxies on the other hand display weaker or 
even absent central point sources and only narrow emission 
lines are viewed through the dusty torus that only allows 
access to the NLR and obscures emission from the BLR and 
accretion disk. However, assuming orientation is the only 
difference means Sy Is and Sy 2s should display virtually 
the same host galaxy properties since the torus only affects 
the very central regions and not the galaxy-wide properties. 
With our focus on the main sequence, this means Sy Is and 
Sy 2s should not separate out in Figure 1. 

Indeed, both Figure 1 and Table 1 suggest this is the 
case. Sy Is compared to Sy 2s in each region show a percent 
difference of +1, +10, -9, -3, and -0 per cent. Assuming 
Poisson statistics, these differences are all well within the 
la error bars. Only the 10 and 9 per cent differences for 
the inside the MS and l-2a regions have a marginal la 
significance. 

To investigate further, we compared the A log SFR dis¬ 
tribution for both Sy Is and Sy 2s and Figure 4 displays 
their histograms. Apart from the increased absolute num¬ 
bers of Sy 2s, Sy Is and Sy 2s have similar distributions of 
A log SFR. Using a K-S test to test whether they are drawn 
from the same parent population, we find a p-value of 0.36, 
again indicating Sy Is and Sy 2s are similar in terms of their 
A log SFR. This is in agreement with Koss et al. (2010) who 
found no difference in u — r colors between broad and narrow 



Figure 4. KDE of A log SFR for Sy Is (solid, blue line) and Sy 
2s (dashed red line) showing the similarity between the two. 


line AGN. This is in disagreement with previous studies (e.g. 
Heckman et al. 1989; Maiolino et al. 1995; Buchanan et al. 
2006) that found Sy 2s reside in more highly star-forming 
galaxies. However, these samples and conclusions are based 
on smaller samples as well as different selection criteria. 
Samples selected using optical or mid-infrared emission will 
inherently be influenced by the level of star-formation in the 
host galaxy and bias samples towards higher star-forming 
objects, especially for Sy 2’s that are fainter at these wave¬ 
lengths due to obscuration. 

5.3 Correlation of A log SFR with host galarxy and 
AGN properties 

Given the large percentage of AGN host galaxies below the 
main sequence compared with normal, main sequence galax¬ 
ies, we examined the relationship between A log SFR and 
various AGN and host galaxy properties. Because our en¬ 
tire sample consists of AGN, the immediate reaction is to 
assume the AGN has influenced star formation in the host 
galaxy through some mechanism and slowed it down. From 
this scenario, the expectation is for more powerful AGN to 
have a greater effect on the host galaxy and occur further 
from the main sequence. To test this, we binned the sam¬ 
ple according to the regions described in Section 5.1. Within 
each bin we calculated the mean 14-195 keV luminosity (Lx) 
since Winter et al. (2012) found the 14-195 keV luminosity 
to be a very good probe of the overall bolometric luminosity 
and strength of the AGN. 

The top left plot of Figure 5 shows the resulting rela¬ 
tionship between A log SFR and Lx. We found no clear cor¬ 
relation between the strength of the AGN and A log SFR. 
Over the entire range of A log SFR, the mean Lx only 
changes by < 0.4 dex with a large spread in each region. 
We calculated the Spearman rank correlation coefficient (ps) 
and used bootstrap analysis to determine the 95 per cent 
confidence interval. We found ps = —0.1 with a 95 percent 
confidence interval of -0.3-0.1 consistent with ps = 0. This 
would seem to argue against the AGN having any effect 
on star-formation in the host galaxy. However Hickox et al. 
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(2014) argues that AGN variability will smear out any in¬ 
trinsic correlations between star-formation and AGN activ¬ 
ity due to the much shorter timescales of AGN activity com¬ 
pared to star-formation. The ultra-hard X-rays used to cal¬ 
culate the luminosity for the BAT AGN presumedly origi¬ 
nate very near the SMBH and represent an instantaneous 
strength while the SFRs are averaged over ~100 Myr. The 
null correlation and large scatter we see between Lx and 
A log SFR then is most likely a product of the large vari¬ 
ability that AGN typically exhibit. 

The middle plot indicates there is a positive correlation 
between Lcrey/LMia and A log SFR. Lorey/LMiR is the ra¬ 
tio of the luminosity of the greybody component to the lu¬ 
minosity of the MIR power-law component used in our SED 
model. The MIR power law luminosity strongly correlates 
with the AGN luminosity (Shimizu et al. 2015, in prepa¬ 
ration) and is assumed to be produced by the AGN, while 
the greybody luminosity is assumed to be a product of star 
formation. The ratio of their luminosities is a measure of 
which component dominates the SED. Far below the MS, 
the IR SED should be dominated by the AGN, while the 
IR SEDs of MS galaxies as well as those above the MS are 
dominated by star formation. We expect this trend due to 
the null correlation of Lx with A log SFR. Since Lmir is 
strongly correlated with Lx, a decrease in A log SFR is then 
mostly due to a decrease in SFR which was directly calcu¬ 
lated from Lgrey. We again calculated ps finding ps = 0.6 
[0.46-0.7] where the range in brackets is the 95 per cent 
confidence interval determined using a bootstrap analysis. 
This strongly suggests a real positive correlation between 
Lcrey/LMiR and A log SFR and confirms that far below the 
main sequence the IR SED is most likely dominated by the 
AGN. 

The greybody dust temperature also suggests a positive 
correlation with A log SFR with galaxies at larger A log SFR 
having a higher temperature (Figure 5, right). This has been 
observed before (e.g. Magnelli et al. 2014) and can easily 
be explained given that an increase in SFR increases the 
number of OB stars that produce the UV photons to heat 
the dust. We found ps = 0.6 [0.45-0.7], very similar to the 
correlation with Lcrey/LMiR showing that the true property 
determining an AGN host galaxy’s location within the main 
sequence diagram is star-formation rather than the strength 
of the AGN. 

Both of these effects can also be seen in Figure 6 where 
we plot the mean SEDs for each region after normalizing 
to the 12 qm flux density. Because some of the observed 
SEDs contain upper limits especially at the longest wave¬ 
lengths, we use the Kaplan-Meier product-limit estimator 
(Feigelson & Nelson 1985), a maximum likelihood estimate 
of the distribution function, to calculate the mean and its 
standard error. Using only the detected flux densities would 
bias the mean towards larger values, with an increasing bias 
at longer wavelengths as the number of upper limits in¬ 
creases. This would then lead to an artificial flattening of 
the SED. 

The peak of the SED is seen to shift to longer wave¬ 
lengths as sSFR decreases, indicating colder temperatures 
(i.e. Wien’s Displacement Law) while the overall ampli¬ 
tude decreases along with a flattening of the slope of the 
SED between 4.6 and 160 qm. In their analysis of Spitzer 
IRS spectra for Seyfert galaxies, both Wu et al. (2009) and 



Figure 6. Mean SED of the sources in each region of the main 
sequence plot. The individual SEDs were first normalized to the 
12 |xm flux density. Error bars represent the standard error of the 
mean normalized flux. 

(Baum et al. 2010) found that AGN host galaxies with lower 
amounts of star formation display bluer SEDs in 15-40 pm 
regime, in agreement with our mean SEDs and an increase in 
the AGN contribution. We also see a general increase in the 
W1/W2 WISE color moving towards lower sSFR, an indica¬ 
tion that the host galaxy (in particular older stars) begins to 
contribute to the MIR, similar to what is observed in many 
quenched galaxies. The long wavelength IR regime, though, 
seems to be completely unaffected by the sSFR. All five of 
the SEDs essentially display the same slope and relative flux 
density at 250, 350, and 500 pm. Indeed a comparison of the 
BAT AGN SPIRE colors with the HRS colors shows there 
is little difference between the two samples (Shimizu et al 
2015, in preparation), verifying the same process (i.e. star 
formation) is producing most of the long wavelength emis¬ 
sion. 

5.4 Host Galaxy Morphology 

Koss et al. (2011) closely analyzed the host galaxy mor¬ 
phologies of the BAT AGN, finding that at all stellar masses 
a larger percentage of AGN are hosted by spiral galaxies 
compared to a matched sample of normal galaxies. Using 
these morphologies, we can assess if we observe a change 
in the host galaxy morphology as a function of sSFR. 
Koss et al. (2011) classified the BAT AGN into three cat¬ 
egories: spirals, ellipticals, and intermediate based on the 
results from the Galaxy Zoo Project (Lintott et al. 2008). 
Each galaxy was independently classified numerous times by 
the public. A spiral or elliptical morphology was chosen for 
the galaxy if > 80 per cent people selected the type, or else 
intermediate was chosen. Mergers were defined in the same 
way as Patton & Atfield (2008) and Koss et al. (2010), re¬ 
quiring a projected distance of at most 30 kpc and a radial 
velocity difference of < 500 km s“^ between the galaxy and 
its companion. We used the same method for the GOLD 
GASS sample as well for the spiral/elliptical/intermediate 
classification. Merger classifications for the GOLD GASS 
sample were determined in Saintonge et al. (2012) where 
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Figure 5. Correlations between different properties of the galaxies/AGN as a function of AlogSFR. left: 14—195 keV luminosity, middle: 
ratio of the MIR power law luminosity and greybody luminosity, right: Dust temperature. For each property, we binned the sources 
according to whether they were above, inside, or below the main sequence (split into 3 separate regions). For each MS region we 
calculated the mean AGN/galaxy property. Error bars are the standard deviation within each bin. 


they visually classified each galaxy as a merger/interaction 
if it had a nearby companion (< 1') and/or evidence of 
a disturbed morphology, tidal tails, etc (see Appendix A 
of Saintonge et al. (2012)). Each galaxy that looked like a 
merger was given a merger rating from 2-5 with 2 repre¬ 
senting galaxies that only had a nearby companion and a 5 
representing galaxies with very strong signs of a merger. Af¬ 
ter looking through the images of the COLD GASS mergers 
we decided to exclude all of the galaxies with a merger rating 
of 2 because they would not have been classified as a merger 
using the method from Koss et al. (2011). We chose to only 
compare the BAT AGN morphologies with the COLD GASS 
sample given COLD GASS’s better mass completeness and 
larger SFR depth. 

In the “Total” column of Table 2 and Figure 7 we outline 
the total fraction of galaxies that are spirals, intermediates, 
and ellipticals over the whole sample. Spiral galaxies domi¬ 
nate (~ 60 per cent) the population of both the BAT AGN 
and GGS similar to what was found in Koss et al. (2011) 
whereas for the CGA they represent ~ 50 per cent of the 
population and for the GGL they are 35 per cent. There are 
virtually no spiral galaxies in the quiescent group (CGQ). 

On the other hand the GGQ are dominated by ellipti¬ 
cal galaxies (~ 60 per cent). The GGL are ~ 26 per cent 
elliptical, CGA 20 per cent, and the BAT AGN and CGS 
10 per cent elliptical. Overall in terms of the whole popula¬ 
tion the BAT AGN more closely resemble the morphology 
distribution of the GGS rather than the CGA. However a 
closer look at the sSFR for each morphology shows that the 
BAT AGN and CGA are more similar. This is shown in Fig¬ 
ure 7, right panel, where we plot the average sSFR for each 
morphology. BAT AGN spirals and intermediates show de¬ 
creased levels of sSFR compared to the CGS and more in 
line with the CGA. If we examine the distribution of spirals, 
ellipticals, and intermediates as shown in Figure 8, we can 
again see that the BAT AGN and GGA are almost identi¬ 
cal in terms of the fraction of spirals and ellipticals in each 
MS region. The biggest difference seems to occur with in¬ 
termediates; however there is much uncertainty concerning 
the nature of intermediates given they are simply defined as 


objects where there was no consensus on whether it was a 
spiral or elliptical. 

Interestingly for both the BAT AGN and CGA, ellipti¬ 
cals show increased levels of sSFR compared to those within 
the CGS and CGQ groups. This could be an indication 
that AGN within early-type galaxies actually stimulate star- 
formation rather than quench it. With such large error bars, 
though, the significance of this result is unclear. 

In agreement with Koss et al. (2010) we find that the 
merger fraction of the BAT AGN (28 per cent) is much 
higher than that of the GGS (3 per cent) and CGA (3 per 
cent). There are no mergers in the CGL and GGQ samples. 
Mergers have long been known to increase the SFR in the in¬ 
teracting galaxies (Sanders et al. 1988) as evidenced by the 
majority of ultra-luminous infrared galaxies (ULIRGS) in¬ 
volved in one. Elbaz et al. (2011) found that the nearly all of 
the local galaxies above the main sequence are IR-compact 
starbursts that were most likely triggered by a merger. Thus, 
we expect to see a large fraction of mergers occurring above 
the main sequence. Indeed, this is the case especially for the 
GGS where 34 per cent of the mergers are above the MS 
and the rest are inside the MS. Further for both the BAT 
AGN and CGS, mergers have the greatest average sSFR 
(Figure 7). This could also explain the extremely low merger 
fraction for the CGA. If mergers are associated with highly 
star-forming galaxies, then emission line ratios would indi¬ 
cate a star-forming galaxy rather than an AGN as any AGN 
signatures would be overwhelmed. With ultra-hard X-rays 
though, we can peer through the obscuring gas and detect 
the AGN. 

In Figure 9, we plot the fraction of all galaxies of a sam¬ 
ple in a MS region that is a merger. This is different than 
Figure 8 where we plot the fraction of all mergers that is 
in each region. Here we can see that the merger incidence 
rate rises with A log SFR as expected. Above the MS, the 
GGA show the highest merger rate at roughly 80 per cent, 
but this is based on small number statistics and is reflected 
in the large error bars. The BAT AGN have a 50 per cent 
merger rate and GGS are at 30 per cent. For all three sam¬ 
ples the above the MS region represents the largest merger 
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Figure 9. Fraction of galaxies in each MS region that are mergers 
for the different samples.. Errors for the BAT AGN (black stars) 
are calculated assuming Poisson statistics for the merger fraction. 
Errors for the CGA (blue squares) and CGS (red triangles) are 
the standard deviation of the fraction over all 50 representative 
subsamples. 

rate. Inside, and below the MS the BAT AGN clearly show 
much larger merger rates than that of CGS and CGA. In 
fact the merger fraction for the BAT AGN is roughly flat 
from above the MS to the 2~Sa region. This could be evi¬ 
dence that mergers and interactions are important in trig¬ 
gering AGN at low redshift even before the star burst ig¬ 
nites and those below the MS are in fact moving up the 
SFR-M, plane. Another explanation is that these are late 
stage mergers and have moved past the initial burst of star 
formation and are slowly falling off the main sequence. A 
confirmation of this would be if all of the mergers in the 
> Icr region are in the early stages whereas the ones above 
the MS are near coalescence. In either case, the BAT AGN 
mergers occur across nearly the whole SFR-M* plane at a 
high rate compared to local non-active galaxies and optically 
selected AGN. This is in disagreement with recent higher- 
redshift studies of the merger rate in X-ray selected galaxies 
(e.g. Kocevski et al. 2012; Villforth et al. 2014) that find no 
difference in the merger rates between AGN and non-AGN 
galaxies. 


6 DISCUSSION 

6.1 Selection Effects and Model Dependence 

Due to the flux-limited nature of our parent AGN sam¬ 
ple, there is a strong dependence on the X-ray luminosity 
(i.e. AGN luminosity) with distance, limiting the inclusion 
of low luminosity AGN especially at higher redshifts. How¬ 
ever this would only affect our results if they were related to 
either the SFR or M* of the host galaxy. X-ray luminosity 
has been shown to only have a weak correlation with SFR 
(e.g Silverman et al. 2009; Shao et al. 2010; Mullaney et al. 
2012a; Rosario et al. 2012). Furthermore, the correlations 
are positive implying the absence of low AGN luminosity 
in our sample only removes objects with even lower SFR 


that would enhance our result that a large fraction of AGN 
lie below the main sequence. Mullaney et al. (2012a) found 
a very weak correlation of intrinsic AGN luminosity with 
stellar mass where a two order increase in AGN luminosity 
results in only a factor 2 greater stellar mass. This com¬ 
bined with the weak correlation of SFR shows that we are 
not biased towards low sSFR galaxies in our sample. 

The choice of model to fit the SEDs of our galaxies 
could greatly influence the results, given that the SFRs are 
directly calculated from the measured FIR luminosity. It is 
possible that we are assigning too much of the 8-1000 pm 
luminosity to the AGN and underestimating the SFR. We 
have tested this effect in two ways. First, we implemented 
the SED fitting routine, DECOMPIR, described and devel¬ 
oped by Mullaney et al. (2011). DECOMPIR utilizes a set 
of 5 host galaxy templates and an intrinsic AGN template 
to fit the IR SED’s of galaxies. Although there is the option 
to let the parameters of the AGN template vary, given our 
sample, we only used the mean intrinsic AGN SED found 
in Mullaney et al. (2011). Each BAT AGN was fit with each 
host galaxy template allowing the normalizations of both the 
AGN and host galaxy to vary. The host galaxy template that 
resulted in the minimum was chosen as the best fit. SFRs 
were then calculated by integrating the best fit host galaxy 
SED from 8-1000 pm and using Equation 1. Using the DE¬ 
COMPIR SFRs combined with the same stellar masses and 
main sequence relation (defined by Equation 3) we find 6, 
35, 28, 14, and 18 per cent of the AGN in the Above MS, 
Inside MS, 1 — 2 ct Below, 2 — 3 ct Below, and > 3 ct Below 
regions respectively. These fractions are extremely similar 
to those using the Casey (2012) model (Table 1) with only 
a slightly higher incidence of AGN inside the main sequence 
and lower number below it. 

The second way we tested the model dependency was by 
using just the monochromatic 160 pm luminosity as a pure 
SFR indicator. The 60 or 70 pm luminosity has been used ex¬ 
tensively in the literature as a SFR indicator; however, based 
on the modeling of the BAT sample, 70 pm is a poor choice 
due to the AGN contribution especially at higher AGN lu¬ 
minosities (Melendez et al. 2014). 160 pm however seems to 
be relatively free from any AGN contribution. Therefore we 
simply used the 160 pm luminosity to convert to a SFR 
using the relation from (Calzetti et al. 2010) for both the 
HRS sample and BAT AGN and recalculated the main se¬ 
quence based only on the HRS since we do not have 160 
pm data for the HerS sample. This results in 6, 37, 25, 12, 
and 20 per cent of the BAT AGN in the Above MS, Inside 
MS, 1 — 2cr Below, 2 — 3a Below, and > 3a Below regions 
respectively, comparable to both the values in Table 1 and 
using DECOMPIR. This confirms our main result that a 
large fraction of AGN lie below the main sequence and is 
not a product of over-subtracting the AGN component of 
the SED. Three separate techniques for estimating the SFR 
agree that > 50 per cent of AGN have lower sSFRs than 
normal galaxies for the same stellar mass. 

6.2 Comparison with previous studies 

We have substantiated that our results are not due to selec¬ 
tion effects or are model dependent. The BAT AGN, COLD 
GASS AGN, and HerS AGN span a large range in sSFR 
that extends over two orders of magnitude. Many previous 
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Table 2. Host Galaxy Morphology Distribution in SFR-M* Plane 


Sample 

(1) 

Total 

(2) 

Above MS 
(3) 

Inside MS 
(4) 

1 — 2f7 Below 
(5) 

2 — 3cr Below 
(6) 

> 3f7 Below 

(7) 

BAT AGN 

Spirals 

0.64 

0.06 

0.33 

0.26 

0.21 

0.14 

Ellipticals 

0.10 

0.00 

0.18 

0.18 

0.09 

0.55 

Intermediates 

0.17 

0.05 

0.11 

0.42 

0.16 

0.26 

Mergers 

0.28 

0.09 

0.34 

0.28 

0.22 

0.06 

CG All 

Spirals 

0.43 

0.03 

0.52 

0.20 

0.12 

0.13 

Ellipticals 

0.23 

0.00 

0.03 

0.03 

0.08 

0.86 

Intermediates 

0.34 

0.03 

0.28 

0.16 

0.13 

0.41 

Mergers 

0.02 

0.41 

0.40 

0.00 

0.19 

0.00 

CGS 

Spirals 

0.59 

0.04 

0.74 

0.16 

0.04 

0.03 

Ellipticals 

0.11 

0.00 

0.06 

0.00 

0.11 

0.84 

Intermediates 

0.29 

0.06 

0.60 

0.11 

0.07 

0.18 

Mergers 

0.03 

0.34 

0.66 

0.00 

0.00 

0.00 

CGA 

Spirals 

0.48 

0.02 

0.34 

0.25 

0.14 

0.25 

Ellipticals 

0.18 

0.00 

0.08 

0.06 

0.14 

0.72 

Intermediates 

0.35 

0.02 

0.29 

0.28 

0.15 

0.26 

Mergers 

0.03 

0.49 

0.01 

0.00 

0.50 

0.00 

CGL 

Spirals 

0.35 

0.00 

0.21 

0.27 

0.26 

0.27 

Ellipticals 

0.26 

0.00 

0.00 

0.08 

0.17 

0.76 

Intermediates 

0.39 

0.00 

0.00 

0.18 

0.30 

0.53 

Mergers 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

CGQ 

Spirals 

0.03 

0.00 

0.00 

0.00 

1.00 

0.00 

Ellipticals 

0.59 

0.00 

0.00 

0.00 

0.00 

1.00 

Intermediates 

0.38 

0.00 

0.00 

0.00 

0.00 

1.00 

Mergers 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


studies have suffered from low detection rates (< 40 per 
cent) in the FIR and needed to resort to stacking techniques 
(e.g. Mullaney et al. 2012a) to achieve a reasonable dynamic 
range. Our Herschel-BAT AGN sample provides a unique 
opportunity to compare to previous work and whether as¬ 
sumptions made remain valid. In particular we will focus on 
the work from Mullaney et al. (2012a) who used deep Her- 
schel observations of the GOODS-South and GOODS-North 
fields to study the star-forming properties of AGN selected 
from 2 Ms and 4Ms Chandra Deep Field South and North. 
Their AGN sample spans the same luminosity range as the 
BAT AGN, but a much higher redshift range (0.5 < z < 3) 
with an overall detection fraction in the FIR of ~ 40 per 
cent. To enhance their sample to low redshifts, they utilize 
a sample of the BAT AGN, however with much poorer FIR 
photometry from IRAS and again only a detection fraction 
of 40 per cent compared to our 95 per cent and 83 per cent 
detection fraction at 70 and 160 pm respectively, the two 
wavelength bands closest to the longest IRAS bands. 

For the FIR undetected AGN in the GOODS fields, 
they utilized stacking analysis to find the weighted average 
sSFR in three different redshift bins. Using a Monte Carlo 
approach and assuming a lognormal distribution of sSFR, 


they estimated 15 per cent of all (FIR detected and unde¬ 
tected) AGN reside in quenched galaxies, 79 per cent are in 
main-sequence galaxies, and 7 per cent in starburst galax¬ 
ies. While the starburst percentage (i.e. above the main se¬ 
quence) is comparable to the fraction presented here, the 
main sequence and quenched (i.e. below the main sequence) 
percentages seem to be in disagreement. However this is 
merely due to a difference in definition of being inside or be¬ 
low the main sequence. Mullaney et al. (2012a) define star- 
bursts as galaxies that have sSFRs > 3 times that of main 
sequence galaxies and quenched galaxies as ones that have 
10 per cent the sSFR of main sequence galaxies. All others 
are then considered normal star-forming galaxies which in¬ 
clude transitioning galaxies in the “green valley.” Using these 
criteria for our BAT sample, we then find 14 per cent star- 
burst galaxies, 73 per cent main sequence galaxies, and 13 
per cent quenched galaxies. While we seem to see a factor of 
2 greater incidence of starburst galaxies, the fraction of main 
sequence galaxies and quenched galaxies are completely in 
line with Mullaney et al. (2012a). 

To further check the agreement between the two stud¬ 
ies, we simulated samples of AGN at higher redshifts by as¬ 
suming our BAT AGN sample is representative of the AGN 
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Figure 7. Left: The distribution of spirals, intermediates, and ellipticals for all of the galaxy samples. For the COLD GASS samples, 
the ractions represent the mean fraction of all 50 representative subsamples. Icr error bars for the BAT AGN (BA) were calculated 
assuming Poisson statistics while the error bars for the COLD GASS samples represent the standard deviation of the fractions for the 50 
representative subsamples. Right: The average sSFR for the different morphologies within each sample. For the COLD GASS samples we 
randomly chose one representative subsample which for the CGA in this case did not contain any mergers. Error bars are the standard 
deviation of sSFR within each group. 


population at all redshifts. To simplify the process, we also 
assumed every AGN has a host galaxy SED exactly that of 
the SB5 template from DECOMPIR at « = 0. The template 
was then scaled to have the same observed 160 pm flux. For 
each redshift bin probed by Mullaney et al. (2012a) (0.5- 
1.0, 1-2, and 2-3) we simulated 1000 samples of AGN by 
randomly assigning a redshift from that bin to each of the 
AGN in our BAT sample. After a redshift was assigned, 
we simulated the increase in SFR with increasing redshift 
using Equation 13 from (Elbaz et al. 2011). From this, we 
obtained a multiplicative factor to boost the IR SED of the 
AGN for a specific redshift. Finally, after redshifting the 
SED, we employed the same detection wavelengths (100 pm 
for 0.5 < 2 < 1.5 and 160 pm for 1.5 < 2 < 3.0) and 
thresholds (0.8 mjy for 100 pm and 2.4 mjy for 160 pm) to 
determine whether the source would be detected. The de¬ 
tection fractions for all 1000 sets of AGN were averaged to 
find the estimated detection fraction for each redshift bin. 

Overall for the entire redshift range of 0.5-3 we found 
a detection fraction of ~ 35 per cent, very near to the ob¬ 
served one of 42 per cent. However when splitting up into 
the redshift bins we find detection fractions of 71, 30, and 
6 per cent compared to 60, 42, and 29 per cent found by 
Mullaney et al. (2012a) in redshift bins of 0.5-1.0, 1-2, and 
2-3 respectively. This could be a sign of an evolution of 
AGN host galaxies where at higher redshifts, more AGN live 
in main sequence galaxies than what is seen at low redshift. 
The discrepancy could also be caused by the lack of an AGN 
contribution used in our simulations. At 100 and 160 pm, 
there is very little contamination from the AGN, but at the 
rest frame wavelengths being probed (40-67 pm) it could 
have a profound effect. To test this we repeated the simu¬ 
lations, adding in the average AGN SED from DEGOMPIR 
such that 25 per cent of the observed frame 100 or 160 pm is 


due to the AGN. Doing this increases the detection fraction 
to 77, 38, and 10 per cent in the same redshift bins, much 
closer in the middle bin but still far off the 30 per cent in 
the highest redshift bin. However, we have not included any 
AGN luminosity evolution with redshift. Ueda et al. (2014) 
find that Lx is a strong function of redshift especially be¬ 
tween 2 = 0 — 2. This would increase the fluxes especially 
in the higher redshift bins and subsequently increase the 
detection fractions. 


6.3 Comparison with SDSS sample 

All three AGN samples, one X-ray selected (BAT AGN) and 
two optically selected (CGA and HerS AGN) have reduced 
levels of star-formation compared to the main sequence. All 
of these samples, however, are relatively small in number so 
we decided to test a much larger sample of both star-forming 
and AGN galaxies from the SDSS. 

The SDSS DR7 spectroscopic catalog mentioned previ¬ 
ously contains 818,333 unique galaxies with optical spectra 
that were analyzed in a consistent manner by the MPA-JHU 
team (Brinchmann et al. 2004). SFRs and stellar masses 
were measured for every galaxy. Stellar masses are based 
on fits to the hve SDSS photometry using the technique 
described in Salim et al. (2007). SFRs were derived in two 
separate ways. For star-forming galaxies, the SDSS spectra 
were fit using the Chariot & Longhetti (2001) models with 
an additional aperture correction to account for light outside 
the fiber (Salim et al. 2007). For the other classifications, 
emission lines are not reliable either due to low S/N or un¬ 
known AGN contribution (composite and AGN dominated). 
The SFRs for these groups were estimated using the 4000 
Abreak {Dn-fOOO) which has been shown to correlate with 
sSFR, albeit with large scatter (Brinchmann et al. 2004). 
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Figure 8. Comparison of the fraction of spirals, ellipticals, intermediates, and mergers that are in each main sequence region for the 
BAT AGN (BA) and COLD GASS samples (CCS, CGA, GGL, and GGQ). The fractions for the GOLD GASS samples are the average 
fractions over the 50 representative subsamples. Errors for the BAT AGN are calculated assuming Poisson statistics while the errors for 
the GOLD GASS samples are the standard deviation of the fraction for the 50 representative subsamples. 


We restricted the SDSS sample to only those galaxies 
with well measured redshifts {zconf > 0.9) and in the red- 
shift range 0.01 < z < 0.05. The upper limit matches the up¬ 
per redshift of our BAT AGN sample and the lower limit is to 
avoid very nearby galaxies that the SDSS pipeline can shred 
into multiple sources. We further included an absolute mag¬ 
nitude cut of Mz < —19.5 to ensure mass completeness and 
only used galaxies with well measured SFRs (sfr_flag = 0). 
These restrictions resulted in 21,091 star-forming galaxies; 
8,189 low-SFR galaxies; 12,190 composite; and 9,421 AGN- 
dominated systems. We combined the star-forming and low- 
SFR galaxies to form the SDSS normal galaxy sample. Fig¬ 
ure 10 shows the distribution of sSFR for the SDSS normal 
galaxies. There is a clear bi-modality in sSFR that matches 
the bi-modality seen with optical colors and in the COLD 
GASS sample. The threshold between the two populations 
occurs at sSFR ~ 10~^^yr“^ shown by the vertical dashed 
line. We use this threshold to split the SDSS normal galaxy 


sample into an actively star-forming population and passive 
population, hereafter referred to as SDSSsf and SDSSpassive. 

In Figure 11 we plot the number density of the SDSS 
galaxies on the SFR-M* diagram separated into the vari¬ 
ous classifications. In the left panel of Figure 11, we plot the 
SDSSsf galaxies (blue contours) along with the SDSSpassiue 
galaxies. The black solid line is the MS relation found by 
Peng et al. (2010) using the same data while the dashed 
line corresponds to a constant sSFR= 10“^^ yr~^. The 
SDSSsf galaxies follow the Peng et al. (2010) MS relation 
just as the HRS galaxies follow their own MS relation with 
the passive galaxies far below. In the right panel of Fig¬ 
ure 11 we plot the SDSS composite (magenta contours) and 
AGN-dominated galaxies (green contours). The same effect 
is seen as when comparing the BAT AGN, CGA, and HerS 
AGN with the main sequence, CGS, and HerS SF. AGN 
host galaxies systematically have lower rates of star 
formation than normal star-forming galaxies. SDSS 
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Figure 11. SFR—M* diagram for the spectroscopic SDSS DR7 sample, left panel: SDSS non-AGN galaxies split into “SF” (blue contours) 
and “passive” (red contours) based on the sSFR. The black line indicates the sSFR= 10“^^ yr~^ threshold to split the sample into “active” 
and “passive”, right panel: Same as left panel with contours representing the SDSS composite (magenta contours) and AGN-dominated 
galaxies (green contours). Gontours enclose 10—90 per cent of the specified sample in increments of 10 per cent. 



Figure 10. KDE of the sSFR for the SDSS DR7 sample. There 
is a clear bi-modality in sSFR defining the two populations of 
actively star forming and passive galaxies with a transition around 
10“^^ yr“^ (dashed line). 

AGN and composites, BAT AGN, CGA, and HerS AGN 
definitively lie in between the actively star-forming and com¬ 
pletely quenched population signaling these gala:xies are pos¬ 
sibly transitioning from one stage to the next. 

6.4 Implications for Galaxy Evolution and AGN 
Feedback 

In this Paper, we have rigorously shown that Seyfert galax¬ 
ies display lower levels of star formation than that expected 
from the main sequence by comparing normal galaxy sam¬ 


ples with an ultra-hard X-ray selected AGN sample using the 
same methods for measuring the SFR and stellar mass. Opti¬ 
cally selected AGN from the COLD GASS and HerS sample 
also show the same effect even though their SFRs were mea¬ 
sured using a different method. Extending the comparison 
to large numbers with the SDSS further emphasizes the dif¬ 
ference between AGN and non-AGN galaxies and confirms 
our results with the much smaller samples. 

Just because AGN host galaxies have systematically 
lower sSFRs than non-AGN galaxies however does not di¬ 
rectly imply that AGN feedback is taking place. Two sce¬ 
narios can possibly explain our results. 1.) AGN actively 
quench star formation through short outbursts during the 
late life of the galaxy. 2.) AGN are simply the result of the 
availability of cold gas in a galaxy. 

High-mass star-forming galaxies on the main sequence 
are gas rich especially in molecular gas (Saintonge et al. 
2012). This large cold gas supply can fuel both star- 
formation and AGN activity. If we prescribe to the “bathtub” 
model of gas regulation (Lilly et al. 2013), the SFR is simply 
proportional to the mass of the gas reservoir in the galaxy. 
The two processes that regulate the mass are the accretion 
rate from the halo and a wind outflow, which in Lilly et al. 
(2013) is proportional to the SFR. 

In the first scenario we can imagine the AGN signif¬ 
icantly adding to the decrease of gas mass in the galaxy 
through several ways such as halo heating and powerful 
winds. All of these feedback processes work to reduce the 
mass of gas in the galaxy, which then reduces the SFR and 
produces the shift in sSFR we see compared to non-AGN 
star-forming galaxies (Figure 3). However based on Figure 5 
we know the SFR is not connected to the instantaneous 
AGN strength, rather the SFR is correlated with the average 
AGN strength over 100 Myr timescales (Chen et al. 2013). 
So if AGN feedback is working it must be only over relatively 
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short periods of time probably while undergoing a powerful 
outburst. This is supported by findings of molecular out¬ 
flows in powerful AGN (Veilleux et al. 2013; Cicone et al. 
2014) where the mass outflow rate was shown to rise with 
increasing AGN luminosity. It is also supported by the dis¬ 
covery of “voorwerpje” (Lintott et al. 2009; Keel et al. 2012, 
2014), highly ionized clouds at kpc scales around currently 
dormant SMBH that indicate AGN outbursts on timescales 
of ~ 100, 000 years. Tombesi et al. (2015) recently were able 
to show that a fast accretion-disk wind is driving a molecu¬ 
lar outflow in a ULIRG. Within this framework, this means 
all of the BAT AGN lying below the MS have gone through 
at least one powerful phase or possibly more to be able to 
significantly deplete the galaxy of molecular gas. The AGN 
still inside the MS perhaps are still waiting for that outburst 
or haven’t gone through enough to move off the MS. 

In the second scenario, the existence of an AGN is the 
consequence of the large availability of cold gas in high 
mass galaxies. If there is enough cold gas, eventually enough 
will find its way to the centers to trigger an AGN. AGN 
would likely turn on while the galaxy is still on the MS. 
A quenching process unrelated to the AGN, possibly shock 
heating of the accreting halo gas, slowly shuts down star 
formation beginning in the outskirts of the galaxy. This 
would be supported by our findings that the BAT AGN are 
more compact in the FIR than normal star-forming galaxies 
(Mushotzky et al. 2014). The AGN persists as the galaxy 
falls off the MS, eventually turning off as the remaining cold 
gas runs out. This scenario is still supported by the hndings 
of Chen et al. (2013). The long timescale average accretion 
rate is tightly connected to the SFR through the available 
gas reservoir. It would only require that gas accretion onto 
the SMBH not be a smooth and constant process but more 
intermittent and variable which is supported by high reso¬ 
lution simulations (Hopkins & Quataert 2010; Novak et al. 
2011 ). 

There is one component we are currently ignoring how¬ 
ever: mergers and interactions. With 30 per cent of the BAT 
AGN involved In one. In both scenarios mergers could be the 
key Ingredient for funneling cold gas to the nuclei of galaxies 
and igniting an AGN. Schawinski et al. (2014) cites major 
mergers as the mechanism to cause fast quenching in early 
type galaxies whereas late type galaxies are slowly quenched 
through a drying up of their gas reservoir. Perhaps both 
quenching mechanisms are at play in the BAT AGN. The 
majority are contained in massive spirals (Figure 7) where 
AGN feedback (scenario 1) or external processes (scenario 
2) are slowly suppressing star-formation while the ones in¬ 
volved In a merger have had their SFRs and SMBH accre¬ 
tion rates briefly elevated before rapidly falling off the main 
sequence. This could also explain why we still see a large 
fraction of mergers just below the main sequence. 

Which scenario is dominant is still a matter of debate 
and at present we are unable to distinguish between these 
concepts. One thing that is clear though is that massive 
galaxies go through an AGN phase as they fall off the main 
sequence. This is shown in Figure 12 where we plot the frac¬ 
tion of all galaxies that are classified as an AGN/Gomposite 
(blue line) or AGN/Composite/LINER (magenta line, if we 
suppose LINERs also contain an AGN) as a function of dis¬ 
tance from the MS in the GOLD GASS sample. The frac¬ 
tions are the average across all 50 representative subsam- 



Figure 12. Fraction of all COLD GASS galaxies that are AGN 
or Gomposite galaxies (blue star) and AGN, Gomposite, or LIN¬ 
ERs (magenta circles) as a function of distance from the main 
sequence. The fractions shown are the average fractions of all 50 
representative subsamples. 

pies. Both lines peak in the 1 — 2cr and 2 —3a below regions 
with close to 80 per cent of the galaxies in these regions 
containing an AGN or LINER. The question that remains 
is whether these AGN have had a substantial effect on the 
SFR in their host galaxy or they are just “along for the ride” 
off the main sequence. An answer could possibly come from 
surveys using integral field unit (IFU) spectroscopy such as 
MaNGA^^ and GALIFA (Sanchez et al. 2012) that provide 
spatially resolved spectra across the entire galaxy. 


7 SUMMARY AND CONCLUSIONS 

Using consistent measures of both stellar mass and SFR, 
we compared three samples of galaxies, one a set of local, 
non-AGN star-forming galaxies (HRS), one a complete set 
of AGN host galaxies selected from the Swift/BAT cata¬ 
log, and one a mass-selected sample of both star-forming 
and AGN host galaxies (GOLD GASS). From the HRS, we 
constructed our own “main sequence” relation and system¬ 
atically analyzed the location of the samples with respect to 
the main sequence. Our main conclusions are as follows: 

1. AGN host galaxies, both X-ray selected and optically 
selected, systematically lie below the main sequence, indi¬ 
cating reduced levels of star-formation. 

2. After splitting the AGN sample into regions of in¬ 
creased offset from the main sequence of star formation, we 
found no dependence of the offset on hard X-ray luminosity. 

3. Analysing the morphologies of the samples, we hnd 
that while the fraction of BAT AGN that are in spirals most 
closely resembles the star-forming sample, the sSFR’s are 
more closely related to the optically selected AGN. The dis¬ 
tribution of sSFR for spirals and ellipticals in the BAT AGN 
also best match the distribution for spirals and ellipticals for 
optically selected AGN. We find a significant increase in the 

https://www.sdss3.org/future/manga.php 
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merger fraction from 0.1 well below the main sequence up 
to 0.4 near, inside, and above it. The merger fractions for 
the BAT AGN are much higher than those for the COLD 
GASS AGN and star-forming galaxies. 

4. These results are both model independent as well as 
unaffected by our selection criteria. A detailed comparison 
with Mullaney et al. (2012a) does not find any discrepancy 
between the two studies, and evolving our population of 
AGN to higher redshifts agrees well with their detection 
fractions. 

5. Expanding the analysis to the larger SDSS sample of 
galaxies shows the effect of AGN host galaxies lying below 
the main sequence also occurs in larger optically selected 
samples and confirms the previous findings of AGN prefer¬ 
entially occurring in galaxies in transition from star-forming 
to quiescence. 

6. We discussed how the trend in sSFR can be explained 
by AGN feedback that reduces the supply of cold gas in the 
galaxy. This slowly suppresses star-formation through short 
periodic outbursts. It can also be explained if the occurrence 
of AGN is simply the result of the availability of cold gas. 
As star formation is quenched through other processes the 
AGN follows along and eventually fades as the cold gas runs 
out. In either case it is clear AGN are prevalent in massive 
galaxies currently falling off the main sequence. 
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APPENDIX A: EXAMPLE FITS TO THE SED 
OF THE BAT AGN AND HRS 

In this Appendix, we provide to the reader example SED fits 
to both the BAT AGN and HRS as well as a brief description 
of the model and fitting procedure. We chose to fit our SEDs 
with the model described in Casey (2012), which is a combina¬ 
tion of an exponentially cutoff mid-infrared (MIR) power law and 
a single temperature greybody. The greybody has the standard 
form, Fjy oc i/^B[y(T), where /3 is the dust emissivity spectral 
index and Bu{T) is the standard Planck function with temper¬ 
ature T, and has been shown to fit very well the FIR SEDs 
of both normal star-forming galaxies as a whole (Gordon et al. 
2010; Galametz et al. 2012; Auld et al. 2013; Cortese et al. 2014) 
and small star-forming regions within the galaxy (Galametz et al. 
2012; Smith et al. 2012). 

The MIR power law, with the form F\ oc , 

can be thought of as the sum of many hot dust components 
that combine to form an overall power law in the broadband 
SED (Kovacs et al. 2010). Casey (2012) found that this simple 
model fits the observed SEDs of ultra luminous infrared galax¬ 
ies (ULIRGs) very well compared to single greybodies and stan¬ 
dard SED template libraries (e.g. Chary h Elbaz 2001). While 
the physical heating mechanism is different between Seyfert galax¬ 
ies and ULIRGS, the same power law shape is observed in MIR 
spectra of AGN (Schweitzer et al. 2006; Mullaney et al. 2011) and 
motivates the addition of this component for the BAT AGN. The 
equation for the total model is then: 

Fu = + NgreyU^B^iT) (Al) 

For each source we used all wavebands that had at least 
a 5(7 detection. The best fit parameters were found from x^~ 
minimization using the Levenberg-Marquardt algorithm for op- 
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timization. Figure A1 displays four randomly chosen BAT AGN 
SEDs and their la uncertainties (black dots with error bars) along 
with the best fit model (black line). The dashed lines indicate 
the best fit MIR cutoff powerlaw (blue) and greybody (red). As 
shown, the model does quite well for the wide variety of SED 
shapes in the BAT AGN sample. 

Figure A2 is the same as Figure A1 except for four randomly 
chosen HRS sources to show the model works equally well for the 
sources without an AGN. 
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Figure Al. Example SED fits for four randomly chosen BAT AGN used for analysis in this Paper. The observed fluxes are shown 
as black dots with Icr error bars, her upper limits to the flux for undetected wavebands are shown as downward pointing arrows. The 
solid black line is the best fit model given by Equation Al while the dashed blue and red lines indicate the best fit MIR powerlaw and 
greybody components, respectively. 
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Figure A2. Same as Figure A1 except for four randomly chosen HRS sources. 
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